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ABSTRACT 

The  foundation  of  the  theory  of  functionally  graded  plates  with  simply  supported  edges,  under  a 
Friedlander  explosive  air-blast  are  developed  within  the  classical  plate  theory  (CPT).  Within  the 
development  of  the  theory,  the  two  constituent  phases,  ceramic  and  metal,  vary  across  the  wall 
thickness  according  to  a  prescribed  power  law.  The  theory  includes  the  geometrical 
nonlinearities,  the  dynamic  effects,  compressive  tensile  edge  loadings,  the  damping  effects,  and 
thermal  effects.  The  static  and  dynamic  solutions  are  developed  leveraging  the  use  of  a  stress 
potential  with  the  Extended-Galerkin  method  and  the  Runge-Kutta  method.  Validations  with 
simpler  cases  within  the  specialized  literature  are  shown.  The  analaysis  focuses  on  how  to 
alleviate  the  unwanted  effects  of  large  deformations  through  proper  material  selection  and  the 
proper  gradation  of  the  constituent  phases  or  materials. 

Key  Words:  Functionally  Graded;  Dynamic  Response;  Explosive  Blast;  Transient  Response 

1.  INTRODUCTION 

During  combat  situations,  the  structure  of  army  military  vehicles  may  have  to  structurally 
endure  the  effects  of  blast  loading.  Advances  in  functionally  graded  materials(FGMs)  which 
combine  the  properties  of  two  dissimilar  materials  has  been  a  motivating  factor  in  viewing  these 
types  of  materials  as  a  viable  alternative  to  the  current  isotropic  metallic  structures  being  utilized 
in  the  hull  and  armor  plating. 
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FGM’s  are  microscopically  nonhomogeneous  with  thermo-mechanical  properties  which  vary 
smoothly  and  continuously  from  one  surface  to  another.  These  graded  structures  allow  the 
integration  of  dissimilar  materials  like  ceramic  and  metals  that  combine  different  or  even 
incompatible  properties  such  as  hardness  and  toughness. 

In  this  paper,  the  foundation  of  the  nonlinear  theory  of  functionally  graded  plate -type  structures 
under  an  explosive  air-blast  is  developed.  An  approximate  solution  methodology  for  the  intricate 
nonlinear  boundary  value  problem  is  devised,  and  results  that  are  likely  to  contribute  to  a  better 
understanding  of  the  structural  behavior  under  an  explosive  blast  with  beneficial  implications 
towards  their  improved  design  and  exploitation  are  presented. 

2.  BASIC  ASSUMPTIONS  AND  PRELIMINARIES 

The  plate  mid  surface  is  referred  to  a  cartesian  orthogonal  system  of  coordinates  ( x,y,z ),  while  z 
is  the  thickness  coordinate  measured  positive  in  the  upwards  direction  from  the  mid-surface  of 
the  plate  with  h  being  the  uniform  plate  thickness  and  y  is  directed  perpendicular  to  the  x-axis  in 
the  plane  of  the  plate.  See  Fig.  1  below. 

The  nonlinear  elastic  theory  of  FG  Plates  is  developed  using  the  classical  plate  deformation 
Theory  [6].  It  is  also  assumed  that  the  FG  plate  is  made-up  of  ceramic  and  metal  phases  whose 
material  properties  vary  smoothly  and  continuously  across  the  wall  thickness.  By  applying  the 
rule  of  mixtures,  the  material  properties  such  as  Young’s  Modulus,  Density,  and  Poisson’s  Ratio 
are  assumed  to  vary  across  the  wall  thickness  as 

P(z)  =  PcVc(z)  +  P mVm(z) ,  (la) 

In  which  Pc  and  Pm  denote  the  temperature-dependent  material  properties  of  the  ceramic  and 
metallic  phases,  of  the  plate,  respectively  and  may  be  expressed  as  a  function  of  temperature  [7, 
9]  as 
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(lb) 


P  =  P0(PlTA  +1  +  PlT  +  P2T2  +p3t 3) 

Where,  Pq,  P  ] ,  P\,  P3,  and  P3  are  the  coefficients  of  temperature  T(K)  and  are  unique  to  the 
constituent  materials.  Vc(z)and  Vm(z)  are  correspondingly,  the  volume  fractions  of  the  ceramic 
and  metal  respectively,  fulfilling  the  relation 

Vc(z)  +  Vm(z)  =  1.  (3) 

By  virtue  of  (3),  Eq.  (la)  can  be  expressed  as 

P(z, T)  =  [Pc (T)  -  Pm  C T)]VC ( z)  +  Pm(T ) .  (4) 

By  observation,  one  can  deduce  that  for  Vf.(z)  =  0,  P(z,T)  -  Pm(T )  and  for  Vc=l, 
P(z,  T)  =  Pc  (D .  As  a  result,  (z)  e  [0, 1] 

Two  Scenarios  of  the  grading  of  the  two  basic  component  phases,  ceramic  and  metal,  through 
the  wall  thickness  are  considered. 

Case  (1):  The  phases  vary  symmetrically  through  the  wall  thickness,  in  the  sense  of  having  full 
ceramic  at  the  outer  surfaces  of  the  plate  and  tending  toward  full  metal  at  the  mid-surface.  For 
this  case,  Vc  (z)  can  be  expressed  as 

vf(z)  =  f^)Wfl±^)+fz£f('lz£IEk)'l  (5) 

W 2)  2  J  (h/2j  l  2  J 

Where  the  signum  function  is  defined  as 

1,  z  >  0 

sgn(z)  =  <0,  z  =  0  (6) 

-1,  z  <  0 
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N  is  termed  the  volume  fraction  index  which  provides  the  material  variation  profile  through  the 
plate  wall  thickness,  (0  <  TV  <  go)  .  A  pictorial  representation  of  the  distribution  of  the  constituent 
materials  are  shown  below  in  Fig.  2. 

Case  (2):  The  phases  vary  non-symmetrically  through  the  wall  thickness,  and  in  this  case  there  is 
full  ceramic  at  the  outer  surface  of  the  plate  wall  and  full  metal  at  its  inner  surface.  For  this  case, 
Vc(z)  can  be  expressed  as 


VAz)  = 


'h  +  2z\N 
v  2  h  , 


(7) 


Below  is  a  pictorial  representation  of  the  antisymmetric  case  shown  in  Fig.  3. 

It  should  be  noted  that  in  contrast  to  case  (2),  where  there  exists  coupling  between  stretching 
and  bending,  such  coupling  is  not  present  for  the  symmetric  case  (1).  Also,  for  the  purposes  of 
simplicity  the  Poisson’s  ratio  will  be  assumed  to  be  constant  throughout  the  plate  structure.  From 
Eqs.  (l)-(7),  the  effective  material  properties  of  a  FG  plate  can  be  expressed  for  the 
asymmetric  case  as 


[E(z,  T),p(z,  T ),  a(z,  T)]  =  [Ecm  (T),  pan  ( T ),  acm  (T)]| 

[Em(T),pm  ( T),am  (T)] 


^  h  +  2  z^N 


2  h 


(8) 


And  for  the  symmetric  case  as 


[E(z,  T ),  p(z,T),a(z,T)]  =  [Ecm  (T),  pcm  (T),  acm  (T)] 

[Em(T),pm(T),am(T)] 


f  _  \N c *  ,  f  _ \N 


\hllj 


1  +  sgn(z) 


-  z 

vV2  j 


^l-sgn(z)N 

l  2  ) 


(9) 


v(z,T)  =  v(T) 


(10) 
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Where, 


Ecm  Ec 


Em  >  Pcm  Pc  Pm 


(11) 


3.  KINEMATIC  EQUATIONS 

3.1  The  3-D  Displacement  Field 

Consistent  with  the  classical  plate  theory  [6],  the  distribution  of  the  3-D  displacement 
quantities  through  the  wall  thickness  can  be  expressed  as 

u  =  u0-zw0x,  v^v0-zw0y,  w  —  vv0  (12a-c) 

Within  these  equations,  (, u,v,w )  are  the  3-D  displacement  quantities  along  the  ( x,y,z ) 
directions,  respectively.  While,  Uq,  Vq,and  up  are  the  2-D  displacement  quantities  of  the  points 
on  the  mid-surface. 

3.2  Non-Linear  Strain-Displacement  Relationships 

The  nonlinear  strain  displacement  relationships  across  the  plate  thickness  at  a  distance  from  the 
mid-surface  take  the  form  [2,  3,  6] 

£ xx  =  u,x  w,x  ^2,  ^ yy  ~  v,y  w,y  ^£xy  —  Y xy  ~  u,y  ~^^,x  w,x^,y  >  (13) 


£xz  syz  H 

Substitution  of  Equations  (12a-c)  into  Equations  (13)  results  in  the  strain  measures  across  the 
plate  thickness  in  terms  of  the  2-D  displacement  quantities  of  the  mid-surface  of  the  plate 
expressed  as 


£xx 

iy°)] 

°xx 

°xx 

Gyy 

>  =  < 

sW) 

byy 

>  +  z< 

b  yy 

J  xy  _ 

v(0) 

I'xy  J 

v(1) 

V  *y  J 

(14a) 
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Where, 


{^S}  *§?  r%)}  =  {u,x+wlx12  v,y+w2y/ 2  Uty+VtX+WtXWty)  (14b) 

i^xx  ^yy  /xy}~i~Wi 0,xx  ~wO,yy  ~2.Wq  xy}  (14c) 

In  the  above  expressions,  ,  £% ,  y% ) ,  are  referred  to  as  the  membrane  strains  and 

,  £yy ,  /-y, )  are  referred  to  as  the  flexural  bending  strains  which  are  also  known  as  the 
curvatures. 

4.  CONSTITUTIVE  EQUATIONS 

The  stress-strain  relationships  for  a  state  of  plane  stress  is  expressed  as  [10] 


®  XX 

'On 

2l2 

0  ' 

Sxx  -oc{z,T)AT 

°yy 

>  = 

Q\2 

Qll 

0 

'  Syy  -a(z,T)AT 

axy . 

0 

0 

266  _ 

/xy 

cr xz  avz  ® zz  ^  2^) 


The  material  stiffnesses,  2,-,  (z),(/  =  1,2,6)  are  given  by  [9,  10] 


<a6,e26)  =  o  (,«) 


1  -v*(T) 


1  -v*(T) 


2[1  +  v(T)\ 


The  standard  force  and  moment  resultants  of  a  plate  are  defined  as 

•/?/  2 


rn/z 

(Ny  ,M jj)  =  J  ( Jij (1,  z)dz  ,  ( i,  j  =  x,y,xy) 


(17) 
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With  the  use  of  Equations  (13) — (17),  the  stress  resultants  and  stress  couples  are  related  to  the 
strains  by  [3] 


\n}-{Nt] 

'[A]  [B] 

v(°>i 

{M}-{M]T\ 

*  — 

[ B ]  [D] 

Wa'ij 

(18) 


In  which  [A],  [B],  and  [D]  are  the  respective  in-surface,  bending-stretching  coupling,  and 
bending  stiffnesses.  For  the  case  of  symmetric  FG  Plates,  [B]=0,  since  there  is  no  bending¬ 
stretching  coupling.  The  global  stiffness  quantities,  A;-y ,  Brj ,  and  D;/- ,  (i,j=l, 2,6)  are  defined  as 


( Ay ,  By ,  D,J )  -  Qij  (1,  z,  z  z  )dz  ,  (i,  j  =  1,2,6) 


h/2 


(19) 


In  view  of  Equations  (8),  (9),  (10),  (16),  and  (19),  the  global  stiffness  quantities  are 
determined  as 


[(An,  ^22  )>  (#1 1  ’  B22 ),  (D\  i ,  Dyj )] 


-2— (£„&>,&,), 
l-v2 


(A\2 ,Bi2 , £>12 )  - - (Ei ,E2,E3),  (20a-c) 

l-v2 

(A66’B66’D66)  =  — - -(EbE2’El) 

2(1  +  v) 

where  for  the  antisymmetric  case, 

El=Ecnih\l  +  \l(N  +  m,  E2  =  Ecmh2 [l/(/V  +  2)-l  KIN  +  2)\  (21) 

E3  =  Emh 3  / 1 2  +  Ean  h3  [1  /( N  +  3)  - 1  /( N  +  2)  + 1  KAN  +  4)] 

And  for  the  symmetric  case, 

Ex  =Eanh[l  +  l/(N  +  l)],  E2  =0,  E3  =  Ecmh3[l/12  +  1/ 4(A  +  3)]  (22) 
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5.  GOVERNING  EQUATIONS 


Hamilton’s  principle  is  used  to  derive  the  equations  of  motion  and  the  boundary  conditions.  It 
is  formulated  as 


SJ  =  s[l(U+V-K)dt  =  0,  (23) 

Jr0 

Where  t0  and  tq are  two  arbitrary  instants  of  time.  U  denotes  the  strain  energy,  V  denotes  the 
work  done  by  surface  tractions,  edge  loads,  and  body  forces,  and  K  denotes  the  kinetic  energy  of 
the  3-D  body  of  the  structure,  while  8  is  the  variational  operator.  In  Equation  (23),  the  strain 
energy  is  given  by 


SU 


rr.  vv  S'  xy  )dzd&0 


(24) 


where  Qq  denotes  the  mid-surface  area  of  the  panel.  The  work  done  by  external  loads  is 
expressed  as 


SV  =  j*  -  Pt(x,  y)Sw(x,  y,  h  /  2) SCIq  -  |*  I* 

jQn  '  '  J  xJ- 


h 

a*yS^’dzdx- 

— 


u 


y  2 


'j  (j^Sudzdy 


(25) 


In  the  above  expression  Pt  (x,  y  )  is  the  distributed  force  at  the  top  surface  (z  =  h/2) ,  (cr^. ,  f7  vv ) 


are  the  specified  stress  components  along  the  plate  edges,  and  (Su,  Sv)  are  the  virtual 
displacements  along  the  normal  and  tangential  directions,  respectively,  along  the  plate  edges. 
Considering  only  the  transversal  inertia  of  the  structure,  the  variation  in  kinetic  energy  is 
expressed  as 

SK  =  I*  iQwSwdClQ  (26) 

Jq0 


Where, 
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(27) 


h  = 


h 

f\p(z)dz 


2 


where  p(z)  is  the  mass  per  unit  volume. 

Considering  Equation  (23)  in  conjunction  with  Equations  (24)-(27) ,  along  with  the  constitutive 
equations  (15),  the  strain-displacement  relationships,  equations  (14),  and  carrying  out  the 
integration  through  the  thickness,  integrating  by  parts  whenever  feasible,  using  the  expression  of 
global  stress  resultants  and  stress  couples,  while  retaining  only  the  transversal  load,  transverse 
inertia,  and  transverse  damping  results  in 


N xy,y)^iQ  +  (  ^yy,y 


N xy,x^^’0  ~Sr[M xx,xx  +  2*^7 xy,xy  +  ^yy,yy 


xx  w0 ,  x  +  N  xy  w0,  y),x 


(N yyw0,y  +  N xyw0,x),y  +  J0^0  ~  Pt(0  ~  Cw0]Sw0}dQ0 


I  {(^xc  ^ xx )^0  M xx<5wo^x  +  xx,x  xy^yy  +  N xxw0,x 

J  v 


Nx)’w0,y)S['v()\  dy  + 


+  j*  \N  yXSllQ  +  (N  yy  N  yy  )  Si’  Q  M  yySWQ^y  +  (M  yy  ^  y  +  2M  ^y^  +  N  xy  WQ  ^ 


Nyy  VV’(j5  y  )  <5tV0  ]  J.V  'jdt  =  0 


(28) 


Invoking  the  arbitrary  and  independent  character  of  variations  du$,  Svq,  Swq  ,  <5w’q  x ,  and  £>vvq  v 

one  obtains  the  equations  of  motion  and  as  a  by-product  the  boundary  terms  or  conditions.  This 
results  in  three  equations  of  motion  in  terms  of  the  global  stress  resultants  and  stress  couples  and 
four  boundary  conditions  along  the  plate  edges.  These  equations  of  motion  and  boundary 
conditions  can  be  respectively  expressed  as 
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^0  :  Nxk,x  +  Nxy,y  -Q 

(31a) 

&0  •  Nyy,y  +  ^ xy,X  ~  ® 

(31b) 

M  xx,xx  +  xy,xy  +^yy,yy  +  ^xxw0,x 

OWq 

Nxyw0,X)y+Pt-Cw0=I0W0 

+  N xy  WQ ,  y  )  x  +  yy  w0,  y  + 

(31c) 

The  boundary  conditions  become 

Along  the  edges  x  =  0,  Lj 

o 

ll 

B 

or 

Uq  =  0 

(35) 

o 

ll 

or 

o 

ll 

o 

- 

(36) 

o 

II 

B 

or 

o 

II 

H 

X 

(37) 

M XX, X  +  xy,  yy  +  ^  XX  w0,x  +  ^  xy  w0,  y  0 

or 

J 

II 

o 

(38) 

Along  the  edges  y  =  0,  L2 

o 

ll 

& 

5: 

or 

Uq  =  0 

(35) 

o 

ii 

£ 

or 

O 

II 

o 

- 

(36) 

o 

II 

S: 

5 

or 

L 

II 

o 

(37) 

M yy,- y  +2M xy,xx  ^xyw0,x  +  ^  yyw0,y  ® 

or 

J 

II 

o 

(38) 

For  the  case  of  all  edges  simply  supported  and  freely  movable  the  boundary  conditions 

are  as 

follows: 

H’O  =  M  ^  =  N ^  =  0,  Nxx  =  N*xx  on  x  = 

o' 

(40) 

w0  =Myy  =  Nyx  =  0,  Nyy  =  N*^  on  y  =  0,  L2 
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It  should  be  mentioned  for  clarification  sake  that  for  compressive  edge  loading, 

AT*  =  -AT0  and  AT*  =  -AT0 
iy  xx  1  y  xx  ’  1  y  yy  iy  yy  • 

AIR-BLAST  LOADING 

With  the  ever  increasing  demands  for  increased  safety  for  the  warfighter  in  the  field  to  operate 
structurally  sound  vehicles  in  the  event  of  an  IED  or  some  other  type  of  explosive,  it  is 
imperative  that  an  understanding  of  the  structural  response  of  various  components  within 
military  combat  vehicles  under  an  explosive  blast  be  understood  so  that  measures  can  be  taken 
from  a  design  standpoint  to  ensure  the  durability  and  survivability  of  these  components.  To  begin 
to  achieve  this  understanding,  the  type  of  explosive  loading  considered  here  is  a  free  in-air 
spherical  air  burst.  Such  an  explosion  creates  a  spherical  shock  wave  which  travels  radially 
outward  in  all  directions  with  diminishing  velocity.  The  form  of  the  incident  blast  wave  from  a 
spherical  charge  is  shown  in  Fig.  4.  Where  Pso  is  the  peak  overpressure  above  ambient  pressure, 
PQ  is  the  ambient  pressure,  ta  is  the  time  of  arrival,  tp  is  the  positive  phase  duration  of  the  blast 
wave,  and  t  is  the  time.  The  waveform  shown  in  Fig.  4.  is  given  by  an  expression  known  as  the 
Friedlander  equation  and  is  given  as 

pt  (0  =  (pso~  po )[!  -  0  -  ta )  /  tp  ]  exp  [-ait  -ta)/  tp  ]  (41) 

Where, 

Pso  =1772/Z3 -114/Z2  +108/Z  (42) 

1/3 

In  equations  (41)  and  (42),  Z  is  known  as  the  scaled  distance  given  by  z  —  R/W  with  R 
being  the  standoff  distance  in  meters  and  W  being  the  equivalent  charge  weight  of  TNT  in  terms 
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of  kilograms.  Also,  a  is  known  as  the  decay  parameter  which  is  determined  by  adjustment  to  a 
pressure  curve  from  a  blast  test. 

For  the  conditions  of  standard  temperature  and  pressure  (STP)  at  sea  level,  the  time  of  arrival 
ta  and  the  positive  phase  duration  t p  can  be  determined  from  [4] 

t/tx  =R/Rl  =(W/Wl)1'3  (43) 

Where  t j  represents  either  the  arrival  time  or  positive  phase  duration  for  a  reference  explosion  of 
charge  weight  W\ ,  and  t  represents  either  the  arrival  time  or  positive  phase  duration  for  any 
explosion  of  charge  weight  W.  The  determination  of  the  standoff  distance  for  any  charge  weight 
W  follows  a  similar  reasoning.  The  application  of  these  relationships  is  known  as  cube  root 
scaling.  It  should  be  understood  that  in  applying  these  relationships  that  the  standoff  distances 
are  themselves  scaled  according  to  the  cube  root  law. 

6.  SOLUTION  METHODOLOGY 

To  satisfy  the  first  two  equations  of  motion,  equations  (31a,b),  a  stress  potential  will  be  utilized 
which  allows  the  in-plane  stress  resultants  to  be  expressed  by  letting 

^ xx  ~  ^,yy  ’  Ny  —  fi.xx  >  ^  xy  ~  ~<P,xy  ■  (44) 

The  third  equation  of  motion,  equation  (31c),  can  be  expressed  in  terms  of  two  unknown 
variables,  the  stress  potential  <p  and  the  transverse  displacement  w0.  To  accomplish  this,  a  partial 
inversion,  of  equation  (18),  the  details  of  which  are  not  presented  here,  needs  to  be  carried  out. 
Performing  a  partial  inversion  results  in  [4] 


id"))  ! 

[A*]  [5*] 

(Wl-{iVr) 

-[B*f  [£>*] 

>))  J 

(45) 
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Where, 


[A*]  =  [A]-\  [5*]  =  — [A]-1[B]  ,  ~[B*f  =[B][A]~l 

[D  ]  =  [D\  -[B] [A]”1  [B]  (46) 


Using  Equations  (44),  (45),  and  (46)  and  simplifying,  Equation  (31c)  takes  the  form 

DV4vt’o  -  (0tXXWQtyy  -  2t/>xyw0xy  +  ^yyWQtXX )  + 1 qWq  +  Cw0  =  Pt  +  (EJ  ~  NT  +  V~M7  (47a) 

Where, 


E1E3  E2  n_  E2 

Ex(\-vZ)  E\ 


(47b, c) 


This  gives  one  governing  equation  with  two  unknowns,  w0  and  0.  One  more  equation  is  needed 
in  terms  of  the  same  unknowns  which  will  give  two  equations  in  terms  of  two  unknowns  which 
can  then  be  solved.  This  will  come  from  the  compatibility  equation.  By  eliminating  the  in-plane 
displacements  from  the  strain-displacement  relationships,  equations  (14)  the  relationship 
between  the  in-plane  strains  and  the  transversal  deflection  known  as  the  compatibility  equation 
can  be  shown  to  be  given  by 

^ xx, yy  +£yy,xx  ~~Yxy,xy  =  W0,xy  ~  w0,xx  wQ,yy  (^8) 

In  view  of  equations  (44)  ,  (45),  and  (46),  the  compatibility  equation  is  expressed  as 

V4</>  =  El(w0'Xy-w0'XXw0,yy)-a-v)V2NT  (49) 


2  2/2  2/2 

In  equations  (47a)  and  (49),  V  =  8  /  dx  +d  fdy  where  V  is  referred  to  as  the  Laplacian 
operator. 

Equations  (47a)  and  (49)  are  the  basic  governing  equations  used  to  investigate  the  structural 

response  of  FG  plates  under  external  excitation  loading.  For  the  purposes  of  this  paper,  from  this 

point  forward,  the  thermal  terms  will  be  discarded.  To  this  end,  to  solve  equations  (47a)  and  (49), 
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the  approach  adopted  from  [2]  will  be  utilized.  In  this  respect,  the  following  functional  forms  are 
assumed  for  Wq  and 


w0(x,  y,  t)  =  wmn  (t)  sin  Amx  sin  juny 


&  ^mn  0 )  CO  ®  ^  mn  0 )  CO ®  ^Pti  3;  C mn  0 )  COS  COS  ^~Pn  3! 

D,Mn  (f )sin  Amxsin  //„  v  +  ^  A*r  y2  x2 


(50a,b) 


Where  /lm  =  nur/  L\ ,  ///;  =  utt/ I^2,m,n  =  1,2,...  are  the  number  of  half  waves  in  the  x  and  y 
directions,  respectively,  and  wmn(t)  is  the  amplitude  of  deflection.  Also, 
Amn  {t),  Bmn  (t),  Cmn  {t),  and  Dmn  it)  are  coefficients  to  be  determined.  By  substituting  equations 
(50a,b)  into  equation  (49),  the  coefficients  Anv}(t),  Bmn(t),  Cnm(t),  and  Dnm(t)  are  determined 
as 


Amn  (0  =  E\  wmn  0)A«  /32A2  ,  Bmn  (?)  =  E\  w^m  (/)/L2  /  32// 2 ,  Cmnit)  =  Dmnit)  =  0  (51) 
As  a  result  of  the  assumed  form  for  the  Airy’s  potential  function,  tf>  the  following  relationships 
expressed  mathematically  as 


rL  2 

* 

rL  1  I  ,  * 

Jo 

,-0  ,Lldy  =  N"L2' 

1  6  Yr  dx  =  N ....  L\ 

J()  r,xx\y =0Ll  yy  1 

(52) 

implies  that  N\  i ,  AL>2  acquire  the  meaning  of  average  in-plane  compressive  edge  loads. 
Moreover,  by  assuming 

Pt  it)  =  Pmn  it)  sin  sin  Amxs\n  jun  y  (53a) 

which  implies  through  integration  of  both  sides  over  the  plate  area  that 

4  f^2  cL\ 

Pmn(t)  = -  Pt  it)  sin  Amx  sin /unydxdy  (53b) 

/.|/.2  J()  Jo 

or  through  integration  gives 
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( m,n )  =  (1, 1) 


(53c) 


Pnm(0  =  16P,  (t)l  /t2, 

and  introduction  of  equations  (50a,b)  and  (51)  into  equation  (47a)  and  retaining  the  resulting 
equation  along  with  the  unsatisfied  boundary  conditions  in  the  energy  functional  and  applying 
the  Extended  Galerkin  technique  results  in  the  following  nonlinear  differential  equation 
governing  the  structural  response  of  FG  plates,  under  external  excitation. 

^’mn  (0  (Onm  Wmn  (t )  +  (Omn  Wmn  (0  ^ mn  (t) (0  —  ^nm  (0  (54) 


Where,  wmn{t )  represents  the  amplitude  of  deflection  of  the  plate  as  a  function  of  time, 
Pmn  (0  =  l6Pt(t)/I07r2  ,comn  =  ^Kmn  //0  is  the  natural  frequency  of  the  FG  plate,  and 

Anm  =  C / 21  QOJmn  is  the  non-dimensional  damping  factor,  and  0.nm  =  E\  (A^n  +  ju%)/\6  .  It 
should  be  noted  that  at  the  center  of  the  plate  (x,  y)  =  (Lj  /2,  L2  /2) ,  (?)  is  equal  to  the 

maximum  deflection  of  the  plate.  In  these  latter  expressions, 


Kmn  =  (£|£3  -  El)KA  (m4  +  2m2 n2 i//2  +  n4V/4)+Nja (m2 +n20>y/2)  (55) 


(1  -v2)Exl\ 


L 


Where, 


*  * 

(//  =  F|  /  L2  is  referred  to  as  the  aspect  ratio  and  ®  =  N  yy  /  is  referred  to  as  the 


compressive/tensile  edge  load  ratio. 

Equation  (54)  is  a  nonlinear  equation  in  terms  of  the  plate  deflections  as  a  function  of  time.  It  is 
interesting  to  note  that  equation  (54)  is  very  similiar  to  Duffing’s  Equation.  To  obtain  the  plate 
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deflections  as  a  function  of  time,  equation  (54)  is  solved  using  the  Fourth-Order  Runge-Kutta 
Method  with  zero  initial  conditions. 

7.  RESULTS  AND  DISCUSSION 

To  validate  the  present  theory,  comparisons  are  made  with  Akay  [1]  who  considered  a  step 
loading  excitation  of  a  simply  supported  elastic  plate  based  on  a  Von-Karman  nonlinear  theory. 
To  make  this  comparison  the  external  excitation,  Pt  (, t ) ,  expressed  as 

Pf(t)  =  P[H(t)-H(t-t0)] 

was  applied.  Hit)  is  referred  to  as  the  Heaviside  Step  function  defined  as 
Hit)  =  \  for  t  >  0  and  Hit)  =  0  for  t  <  0  The  geometrical  and  material  properties  used  were, 

L|  =  2.438  m,  h  =  0.00635  m,  (//  =  1 
Em  =  70.3  GPa,  pm  =  2547  kg  I  m3 ,  vave  =  0.25 

In  addition,  the  constant  pressure  and  time  duration  were  taken  as  P  -  48.82  Pa,  Iq  -  0.2  sec  . 

Since  only  a  metallic  isotropic  plate  was  made  for  comparison,  the  volume  fraction  index  was 

taken  as  N  =  2000  which  implies  fully  metal.  With  this  in  hand,  the  central  deflection  time 

history  is  displayed  in  Fig  5.  The  results  in  this  figure  reveal  close  agreement  between  Akay’s 

finite  element  method  and  the  present  analytical/  approximation  method  employed  here  in  the 

present  analysis. 

To  illustrate  the  present  approach,  a  ceramic-metal  functionally  graded  plate  consisting  of  Ti- 
6A1-4V  and  Aluminum  Oxide  with  the  following  material  properties,  which  were  adopted  from 
[9],  were  considered  for  the  numerical  results  presented. 

Ec  =  320.24  GPa,  pc  =  3750  kg  /  m3 ,  vc  =  0.26 

Em  =  105.7  GPa,  pm  =  4429  kg  /  m3 ,  vm  =  0.2981 

vave  =0.2791 
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The  geometrical  properties  used  for  the  FG  Plate  are  Lj  =lm,  (//  =  L\  / L2  =  1 ,  and  unless 
otherwise  stated  h  =  0.0254  m,  the  decay  factor,  a  =0.5,  and  the  halfwaves,  (. m,n )  =  (1,1).  In 
addition,  the  following  reference  values,  in  Table  1.  were  utilized  to  determine  the  time  of  arrival 
and  positive  phase  duration  [5].  In  Fig  6,  comparisons  of  the  central  deflection  of  the  plate  for 
various  volume  fraction  indexes  is  depicted,  it  can  be  seen  that  the  symmteric  functionally 
graded  case  with  a  volume  fraction  index  of  N=0.5  gives  the  lowest  deflections  as  a  function  of 
time,  as  compared  to  the  fully  metal  isotropic  case. 

Fig  7.  is  the  counterpart  of  Fig  6.  with  the  exception  of  the  effect  of  damping.  The  effect  of 
damping  shows  a  rapid  attentuation  of  the  deflections  over  a  very  short  period  of  time.  It  is 
shown  that  damping  plays  a  major  role  in  the  decrease  of  deflections  of  the  plate. 

In  Fig  8.  it  can  be  seen  that  the  effects  of  various  amounts  of  damping  on  the  central 
deflections  vs  time  attenuates  faster  as  the  amount  of  damping  is  increased  for  a  fixed  volume 
fraction  of  the  constituent  materials.  The  effects  of  the  compressive/tensile  edge  loading  on  a  FG 
plate  for  a  fixed  volume  fraction  index  and  a  fixed  amount  of  damping  on  the  deflection -time 
response  is  depicted  in  Fig.  9.  It  is  clear  that,  for  tensile  edge  loading,  the  magnitudes  of  the 
deflections  decreases  in  contrast  to  compressive  tensile  edge  loading  where  the  deflections  are 
increased. 

Fig  10.  reveals  the  effects  of  the  aspect  ratio  on  the  deflection-time  response  which  indicates 
that  the  deflections  are  larger  for  smaller  aspect  ratios  than  for  larger  aspect  ratios.  From  a  design 
standpoint,  utilizing  larger  aspect  ratio  panels  would  be  more  beneficial  by  give  smaller 
deflections  under  an  explosive-type  excitation,  hence  decreased  stresses. 

Fig  11.  shows  that  as  the  volume  fraction  index  increases  from  a  value  of  zero  (fully  ceramic) 
to  a  value  of  500  (essentially  fully  metal)  that  the  deflections  increse  then  taper  off  as  phases  of 
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the  constituent  materials  approach  fully  metal.  Fig  12.  shows  the  effect  of  the  decay  parameter 
on  the  deflection-time  history  response.  For  small  values  of  the  decay  parameter,  it  can  be  seen 
that  the  negative  phase  of  the  blast  pulse  or  suction  phase  lasts  longer  and  as  the  the  value  of  the 
decay  paremeter  increases  the  negative  phase  endures  less  thus  giving  more  periodic  behavior 
around  the  zero  deflcetion  point. 

8.  CONCLUDING  REMARKS 

A  rigorous  treatment  of  functionally  graded  plates  with  grading  in  the  transverse  direction  has 
been  studied.  Validations  with  a  simpler  transversal  excitation  step  load  found  in  Akay  [1],  who 
took  a  specialized  finite  element  approach,  has  been  made.  It  has  been  shown  that  damping  has 
an  important  effect  when  it  comes  to  the  attenuation  of  the  deflections.  It  has  also  been  shown 
that  other  factors  such  as  the  compressive/tensile  edge  loading,  the  aspect  ratio,  and  the 
symmetry  of  the  transverse  grading  throughout  the  structure  plays  an  important  role  in  the 
deflection-time  history  of  the  structure. 

In  Fig.  6.  it  was  shown  that  functionally  grading  inherently  reduces  the  deflections  when 
compared  to  the  isotropic  metallic  case.  From  a  design  standpoint,  it  would  be  appropriate  at  this 
point  to  state  that  integration  of  functionally  graded  materials  within  plate-type  structures  would 
benefit  the  structural  response  of  the  structure.  Also,  it  should  be  mentioned,  although  not  shown 
here,  that  the  choice  of  the  ceramic  and  metal  constituent  materials  chosen  would  also  have  a 
great  impact  on  the  response  of  the  dynamic  response  of  the  structure. 

The  idea  is  to  reduce  the  stresses  within  the  structure  concerned  here.  By  reducing  the 
magnitude  of  the  deflections,  the  stresses  are  reduced.  It  is  hoped  and  realized  that  this  present 
study  presented  here  will  give  insight  into  some  of  the  factors  that  can  play  an  important  role  in 
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the  structural  response  of  functionally  graded  plates  and  fill  in  some  of  the  fundamental  missing 

gaps  within  this  subject  area. 
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FIGURE  CAPTIONS 


Fig.  1  A  simply  supported  functionally  graded  plate  shown  in  2-D  under  an  explosive  blast. 

Fig.  2  Distribution  of  the  constituent  materials  through  the  plate  thickness  for  the  symmetric 
case. 

Fig.  3  Distribution  of  the  constituent  materials  through  the  plate  thickness  for  the  antisymmetric 
case. 

Fig  4.  Incident  pressure  profile  of  a  blast  wave. 

Fig.  5  Comparisons  of  solutions  of  the  time  history  of  the  central  deflection  under  a  heaviside 
step  load. 

Fig.  6  Implications  of  the  volume  fraction  index  on  the  central  deflection  as  a  function 
of  time,  (i? =0.766  m,  ta  =  0.000408s ,  tp  =  0.001372s ) 

Fig.  7  The  counterpart  of  Fig  6  with  a  damping  measure  of  Aj  ^  =0.1 

Fig.  8  Implications  of  various  amounts  of  damping  on  the  deflection-time  response  of  an 
antisymmetric  FG  Plate.  (N=0.5) 

Fig.  9  The  implications  of  the  compressive/tensile  edge  loading  on  the  deflection-time 
response  of  an  antisymmetric  FG  plate.  ( =0.2,  N=0.5) 

Fig.  10  The  effects  of  various  aspect  ratios  of  an  antisymmetric  FG  plate  on  the 
deflection-time  response.  (N=0.5) 

Fig.  1 1  Implications  of  the  volume  fraction  index  on  the  deflection-time  history 
of  an  antisymmetric  FG  Plate  without  damping. 

Fig.  12  The  implications  of  the  value  of  the  decay  parameter  on  the  time -deflection 
response  of  a  antisymmetric  FG  Plate  without  damping.  (N=0.5) 
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